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. We consider neutrino oscillations in a medium in which scattering pro- 

cesses are blind to the neutrino flavour. We present an analytical derivation 
^ I of the synchronised behaviour obtained in the limit where the average scat- 

tering rate is much larger than the oscillation frequency. We also examine 
■ MSW transitions in these circumstances, and show that a sharpening of the 

. transition can result. 
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I. INTRODUCTION 



^ ■ We examine a system consisting of an oscillating neutrino propagating through a dense 

Q . medium. This falls into the very broad class of problems in which a quantum system cannot 
be considered closed or isolated, due to interaction of the system with an environment or 
Ph! heat bath. In particular we consider two-flavour neutrino oscillations in a medium such that 

p il the system-bath interactions do not distinguish the flavour of the neutrino. 

An example is the density matrix of a single neutrino undergoing — oscillations 
embedded in an electron-positron plasma. In some sense this example is similar to an early 
universe environment, though lacking some features of a realistic physical scenario, such as 
^ \ neutrino pair processes and the inclusion of other neutrino flavours. Another example in 
d I which flavour blind interactions are relevant is the case of — Vr oscillations of neutrinos 
coming from the core of a supernova. In this case, although the indices of refraction for and 
Vj. will be different due to loop graph effects, the collision rates for the neutrinos scattering 
from the medium of nucleons will be the same for the two flavours, to a good approximation. 
However, it is the coherence properties of the system that are most intriguing, and which 
are worth studying independently of the various circumstances in which they may arise. 

An analogue to the neutrino system is that of a double well coupled to a thermal bath, in 
which the neutrino flavour oscillations are replaced by oscillations from the left to the right 
side of the well. Our work provides an analytic treatment of the distinctive behaviour that 
we noted previously in a numerical treatment of double well oscillations [|l|. Since double 
well systems are of general interest in condensed matter physics the type of behaviour 
we examine here may be expected to arise in various situations unrelated to neutrino physics. 

As the medium cannot distinguish the flavour of the neutrino, the only effect of collisions 
with particles in the medium is to change the value of the neutrino's energy (and to entangle 
the neutrino with the environment). This is in contrast with active-sterile oscillations in 
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which colhsions may in some sense be thought of as measuring the flavour of the neutrino. It 
is the consideration of the neutrino energy variable, and the rapid energy changing colhsional 
processes that is crucial, and when the average collision rate is much larger than average 
oscillation rate we flnd that a collective oscillation results That is, for a neutrino in a 
thermally distributed mixture of momentum states, the oscillations are synchronised despite 
being momentum dependent in the absence of collisions. 

We stress that the novel synchronisation effect arises as a result of energy changing 
processes. In particular, the result is distinct from the synchronisation effect of ref. |^ 
which results from non-linear feedback, as we consider here only linear dynamics. 



A. The Quantum Kinetic Equations 

The evolution equations for the neutrino density matrix have been formulated in ref. 
for active-active oscillations, and ref. for active-sterile oscillations. These Quantum 
Kinetic Equations (QKEs) or Bloch equations, take into account both coherent oscillation 
behaviour and potentially decoherent effects of interaction of the neutrino system with an en- 
vironment. The QKEs are essentially Boltzmann equations, generalised to allow for quantal 
coherence between particle species. 

We choose to parametrise the neutrino density matrix as 

p{p) = l[Poip) + P{p)-cr], (1) 

with 

P{p) = P^{p)x + Py{p)y + P,{p)z. (2) 

The QKEs for a single neutrino undergoing elastic scattering in a medium which does 
not distinguish flavour are 

P{k) = J Sk' [r(fc', k)P{k') - V{k, k')P{k)\ + V(A;) X P(fc), 
P^{k) = I d^k'[Tik',k)Poik') ~Tik,k')Poik)], (3) 

where r{k',k) determines the rate at which neutrinos are scattered from momentum state 
k' to state k. Thermal equilibrium of the medium enforces 

r(fc',fc) = e(*='~'=)/^r(A;,A;'). (4) 

The cross product term describes precession of the neutrino flavour, where V(A;) is given by 

V{k) = (3{k)x + X{k)z, (5) 

and the oscillation parameters are deflned as 

/3(A;) = ^sin2^o, 

X{k) = -^cos29o + V,s, (6) 

with 6*0 being the vacuum mixing angle, 6m'^ the mass-squared difference, and Kg the 
effective matter potential resulting from coherent forward scattering. 
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B. Qualitative features for different coupling strength. 



By way of example, we consider the density matrix of a single v^—Pt neutrino propagating 
in an electron-positron plasma. For simplicity, we shall ignore the interaction of the neutrino 
with other u^^t or 17^,- that may be produced in the plasma, and consider only the single 
neutrino and its interaction with electrons and positrons. That is, we consider in isolation 
a particular subset of the various processes that may affect the evolution of the neutrino 
density matrix. 

For our considerations it should suffice to take the following form for the scattering rate, 



T{k, k') 



0.01 



(7) 



as this embodies the correct threshold behaviour, the reciprocity relation (§) and leads to the 
correct total transition rate. We have solved the QKEs @ numerically, and the results are 
presented in the graphs below, for various temperatures. For small temperatures, where the 
average time between collisions is much longer than the oscillation period, we have just the 
superposition of (almost) undamped precession rates at different frequencies. For moderate 
collisions rates we rapidly lose coherence in the neutrino system. This occurs simply because 
the oscillation modes are momentum dependent. Collisional processes transfer quanta across 
the momentum states which quickly get out of phase, and hence result in a damping of the 
amplitude of the oscillations. The most fascinating effect, however, is obtained for larger 
values of F where the average time between collisions is much smaller than the oscillation 
period. In this limit we see a synchronisation of the oscillations of all the momentum states. 
Here, the precession proceeds at at a rate determined by the thermal average of V(/c), and 
damping is eliminated as the rapid collisions force the momentum modes to stay in phase. 
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FIG. 1. The probabiUty that the neutrino is a z^^, where the initial condition is a with a 
thermal distribution of energies. The sohd, dotted, dashed and heavy solid curves correspond to 
T=5, 10,20 and 50MeV respectively, and we have taken 5m'^ = O.OOleV^, sin 2^0 = 1- In each case, 
the time is scaled by the inverse of the thermally averaged oscillation frequency. 
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We define the entropy in the neutrino system as 

d^kTT[p{k) In p{k)] 



(8) 



where the trace is over the neutrino flavour index. The entropy is plotted in flg.2 for various 
colhsion rates. 
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FIG. 2. Entropy plotted for the same parameters as fig.l, where S=s(t)-s(t=0). 

Let us compare this behaviour with the familiar Quantum Zeno behaviour in active- 
sterile neutrino systems. In both cases the rate of entropy growth increases as we increase the 
interaction rate between the neutrino system and the environment, and then decreases again 
if we make the interaction rate much larger than the oscillation rate. In the active-sterile 
case however, this occurs because the "measurement-like" interactions freeze the neutrino 
in its initial state. In contrast, in our case the flavour precession proceeds as though it were 
a single mode of an isolated system. 



II. ANALYTIC SOLUTION FOR LARGE COLLISION RATE 

We now show analytically that a synchronised solution follows from the Eqs. in the 
large F limit. We begin be discretizing the equations, by allowing the momentum to take 
the values fci, /c2, ••••^Af- Expressing the resulting 3N equations in matrix form we have 



V = MV, 



where V is the vector 



— (-Pxl; PxN^ Pyl, PyN , Pzl, Pzn)'^ ■ 



The matrix Ai, has the following block form 



g 


-A 





A 


g 


-B 





B 


g 



(9) 



(10) 



(11) 
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where 



fc5 



A = diag(Ai, An), 

S = diag(/3i, /5n). (12) 

To solve for P(t), we introduce a matrix U which diagonahses M. such that 

V=UQ, (13) 

and 

Q = MdiagQ-W-^WQ, (14) 

where 

A^diag = U-^MU. (15) 
We shall consider time independent A and /3, so that U = 0, and the solution to eg. (p!4D is 

m) = E^^.SiW = Eexp [(A^diag),,t] W.,Q,(0). (16) 

i i 

The procedure may also be adapted for time dependent mixing where A and f3 change 
sufficiently slowly that an adiabatic approximation is valid]^ 

Since wish to find the solution in the large F limit, we first drop the A and (3 terms. Then 
we see that has three zero eigenvalues, with eigenvectors that are linear combinations of 

yi = 4^ (e-^/^, 0, O)"^ , y2 = ^ (O, e-'/^, of and ys = ^ (o, 0, e~'/^Y , (17) 



where e~^^'^ denotes (e^'^^/"^, e"'^^/"^, e"'^^/^). The corresponding eigenvectors of Ai^ are 
xi = [k\ 0, O)^ , X2 = (O, k\ of and ^3 = ^ (o, 0, k'Y , (18) 

with UU^^ = 1 requiring V NN' = J2k k'^G~^^'^ ■ Q All the other eigenvalues are negative and 
proportional to F. To observe this let VL be an eigenvalue of with eigenvector X. The 
eigenvalue may be expressed as 

^ = Er,^^'(^-l)' Vj:X,^0. (19) 



^An adiabatic approximation was applied to the case of partially coherent oscillations in Ref. 0. 
^Note that since M. is neither hermitian nor symmetric, U will neither be unitary nor orthogonal. 
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Since we may choose j such that \Xj/k'j \ > \Xk/kl\, Vfc, it follows that Reft < 0. Thus the 
eigenvectors corresponding to the non-zero eigenvalues will make a negligible contribution 
to the solution in the large F limit (for t > l/F). Then the solution is simply 

Pj(t) ^ e~^^^'^ X const = thermal solution. (20) 

Treating A and 13 as perturbation parameters, we are led to consider the 3x3 matrix 
Vij = xfVyj, where V is the matrix obtained by setting Q = in A^. Using standard 
degenerate perturbation theory, to first order the three eigenvalues become 

±iw, (21) 

where u, the synchronised oscillation frequency, isQ 

= ^/W+W (22) 
and (A) and {(3) are the thermally averaged values 

^ ^ ~ Ek' k'^e-^''^ ' " E.' A:'2e-^VT " ^^^^ 
Note that in the absence of a matter potential, the synchronised precession frequency is just 
uj = {5m? /2k). For time independent mixing the solution becomes, 

Px{k,t) = ne"^/^C2eS2e[A' - Acos{ujt - a)] 
Py{k, t) = -ne~^^'^ As2@ sin(tut - a) 

P,{k, t) = ne~^/^ [A'cIq + Asj^ cos{ut -a)], (24) 

where 

C2e = 

n= ^J^.,IT T.k"P.{k\t = Q). (25) 

That is, we have undamped oscillations with frequency uj at an effective mixing angle of G. 
If the initial condition is such that the neutrino is in the pure flavour state described by 
(1 + cr3)/2 then a = and A' = A = 1. Allowing for general initial conditions, such as an 
initial distribution of phases, we have 

Ek k'PyiK 0) 



tana 
A' 



Efc' fc'2[s2©P,(fc', 0) - C2ePAk', 0)] 
, , S2eEkk^P.{k,0) 



C2eEk' k'^P,{k',0) 

A = 



Ek k'Pyik, 0) Y ^[^ C2e Ek k'P.ik, 0) ^^ ' 



S2e Ek' k'^P^ik', 0)J \ S2e Ek' k'^P,{k', 0)^ 



(26) 



^If one assumes the existence of a synchronised solution, the oscillation frequency must be as given 
in eq.(22), from the self-consistency of eq.(S). 
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III. MSW SHARPENING 



We now consider the case where there is an MSW resonance, but the colhsion rates are 
still flavour blind. For example, imagine that our plasma has a nucleon component, or a 
lepton excess. The effective potential arises at 1 loop level, and is given by M 



V 



3G}m^ 
27r2 



{Np + Nn)\n 



m 



w 



mt 



(27) 



We need only a very small nucleon component for this potential to be signiflcant. For 
example, assuming a mass squared difference of drn^ ~ 10~^eV^, a nucleon density of Nn ~ 
10~'^A'"e would give rise to an MSW resonance at a temperature of roughly 10 MeV. If the 
nucleon density were smaller, the resonance would occur at higher temperature. 

To illustrate the type of effect that could arise, we assume, by way of example, an 
exponentially decreasing effective matter potential. We take for our initial condition a 
neutrino in a deflnite flavour state with a thermal distribution of momentum values, and plot 
below the probability that the neutrino remains in that initial flavour state. In the absence 
of neutrino-environment interaction, each momentum state goes through the resonance at a 
different time resulting in a broad transition region, while for moderate interaction rate the 
oscillations get dephased and the neutrino ends up in an equal mixture of the two flavour 
states. For large interaction rate, we see the effect of synchronised behaviour - all momentum 
states go through the resonance at the same time, resulting in a much sharper transition. 
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FIG. 3. "MSW sharpening". The probabihty that the neutrino remains a z^^, given an initial 
thermal distribution of f^. The dashed, dotted and sohd curves are for T=l, 10 and 50 MeV 
respectively, with 6m? = O.OOleV ^ and sin 2^0 = 0.2. The time is in units of l/[uj(V = 0)]. 

The sharpened MSW resonance when occurs (A) = 0, where the matter potential has 
the value 

_ E.fc2e-^/^c,,„^mV(2fc) 
^efr(res) p,-./T • (28) 
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IV. CONCLUSION 



When rapid scattering processes are flavour blind, a neutrino in a mixture of momentum 
states undergoes synchronised flavour oscillations. This collective, dissipationless behaviour 
has been demonstrated numerically and derived analytically. Under these circumstances 
MSW transitions become sharper as all momentum states go through the resonance together. 
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